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ABSTRACT 

We present Submillimeter Array observations of the Herbig Ae star HD 169142 
in 1.3 millimeter continuum emission and 12 CO J=2-l line emission at ~ 1'.'5 
resolution that reveal a circumstcllar disk. The continuum emission is centered on 
the star position and resolved, and provides a mass estimate of ~ 0.02 M & for the 
disk. The CO images show patterns in position and velocity that are well matched 
by a disk in Keplerian rotation with low inclination to the line-of-sight. We use 
radiative transfer calculations based on a flared, passive disk model to constrain 
the disk parameters by comparison to the spectral line emission. The derived 
disk radius is 235 AU, and the inclination is 13°. The model also necessitates 
modest depletion of the CO molecules, similar to that found in Keplerian disks 
around T Tauri stars. 

Subject headings: stars: individual(HD169142) — stars: pre-main sequence - 
stars: planetary systems: protoplanetary disks 


1. Introduction 

Millimeter-wave observations of disks surrounding pre-main sequence stars are a promis¬ 
ing probe of the planet formation process (Beckwith & Sargent 1996). While most studies 
have concentrated on the circumstellar disks around low mass T Tauri stars, far less is known 
about the disks around their intermediate mass counterparts, the Herbig Ae/Be stars. Better 
characterization of the circumstellar disks around these higher mass stars will help address 
how stellar mass affects the development of planetary systems. 
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Detailed modeling of the spectral energy distributions of Herbig Ae/Be stars has gen¬ 
erally supported the protoplanetary disk paradigm (e.g. Hillenbrand et al. 1992; Chiang et 
al. 2001; Dominik et ah 2003). But out of the large sample of Herbig Ae/Be star disk can¬ 
didates recognized by infrared excess from circumstellar dust (The et ah 1994; Malfait et 
ah 1998), only a handful have been spatially and spectrally resolved through direct imaging 
observations to show clear signatures of rotationally supported disks of Solar System size 
(Mannings & Sargent 1997, 2000; Pietu et ah 2003; Corder et ah 2005). 

The emission line star HD169142 (SAO 186777, MWC 925; spectral type A5) with 
infrared excess (IRAS 18213-2948) was included in the The et ah (1994) catalog of 287 
Herbig Ae/Be candidates. Subsequent photometric and spectroscopic observations filled out 
its spectral energy distribution at infrared and submillimeter wavelengths (Sylvester et ah 
1996; Dunkin et ah 1997; Malfait et ah 1998). At ~ 145 pc distance, HD169142 is notably 
isolated from any molecular cloud or extensive reflection nebulosity, though Kuhn et ah 
(2001) used differential polarimetry to detect a ~ 200 AU (~ T.'5) diameter circular halo of 
scattered light. While some doubt lingers as to the pre-main-sequence nature of HD169142, it 
is among isolated Herbig Ae stars selected by Meeus et ah (2001) and classified as “group lb” 
with a strong mid-infrared component, and a conspicuous absence of 10 /jrn silicate emission 
(Bouwman et ah 2001). This star also has prominent polycyclic aromatic hydrocarbon 
emission features (Habart et ah 2005). 

The spectral energy distribution of HD169142 is well fit by a passive irradiated disk 
model, though the near-infrared flux seems low, perhaps because the inner disk rim is viewed 
polc-on (Dominik et ah 2003), or perhaps because extensive inner disk clearing is underway, 
possibly together with planet formation. The relative isolation of this star from any molec¬ 
ular clouds may indicate an advanced pre-main-sequence age. A circumstellar molecular 
component is detected via relatively strong emission in the CO J=3-2 line, with the narrow 
CO line profile modeled with a near face-on disk (Greaves et ah 2000; Dent et ah 2005a). 

Here we present high resolution observations of HD 169142 from the Submillimeter Ar¬ 
ray 1 (SMA) that clearly resolve a circumstellar disk in Keplerian rotation. We derive the 
basic disk properties from arcsecond resolution images of the 12 CO J=2-l line using radiative 
transfer calculations based on a flared, passive disk model (D’Alessio et ah 2005). 


Whe Submillimeter Array is a joint project between the Smithsonian Astrophysical Observatory and the 
Academia Sinica Institute of Astronomy and Astrophysics, and is funded by the Smithsonian Institution and 
the Academia Sinica. 
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2. Observations 

We used the SMA (Ho et al. 2004) to observe HD169142 on 19 April 2005. Table 1 lists 
the observational parameters. Seven antennas were available in an extended configuration 
that provided a baseline range of 27 to 181 meters and a synthesized beam of 1"6 x l'/O” 
at position angle 20°. The target source was tracked over the hour angle range -1.0 to 3.5. 
Weather conditions were excellent, and system temperatures ranged from 85 to 150 K (DSB). 
The correlator provided 2 GHz of bandwidth in each sideband and was configured to include 
the 12 CO J=2-l line at 230.538 GHz in the upper sideband in a 104 MHz spectral “chunk” 
with channel spacing 0.2 MHz (~0.27 km s -1 ). Complex gain calibration was performed 
using the calibrator sources J1911-201 and J1924-292. Passband calibration was done with 
observations of the strong sources 3c279 and Uranus. The absolute flux scale was set using 
observations of Callisto and is estimated to be accurate to about 10%. All of the calibration 
was done using MIR software, followed by standard imaging and visual analysis using the 
miriad package. 


3. Results and Discussion 
3.1. Continuum 

The lower right panel of Figure la shows the continuum image made from the upper 
sideband visibility data. Fitting an elliptical Gaussian to the visibilities gives a flux of 
169 ± 5 mJy, and a fwhrn size of ~ 079. The fit does not indicate any significant ellipticity. 
The observed flux is consistent with the previous single dish 1.3 mm measurement of 197 ± 
15 mJy (Sylvester et al. 1996) taking account of the uncertainties, which suggests that the 
interferometer is not missing any large scale emission. The 1.3 mm continuum emission arises 
from a compact region, which we associate with the circumstellar disk. The spectral energy 
distribution indicates that thermal dust emission strongly dominates at this wavelength; any 
contribution from ionized gas emission is negligible. 

The dust continuum flux, S„, may be used to estimate the disk mass, M^sk , assuming 
optically thin emission, i.e. M^sk = S U D 2 / k„B u ((T)) where B u is the Planck function. 
We adopt the mass opacity of Beckwith et al. (1990), = 0.02(1.3 mm/A) cm 2 g -1 , 

which includes a gas-to-dust ratio of 100. Assuming an average temperature ( T ) = 30 K, 
appropriate for the bulk of the disk material, the disk mass is 0.02 M©. The factor of few 
uncertainty in k u likely dominates the uncertainty in this mass estimate. This mass estimate 
falls in the middle of the range of estimates for protoplanetary disks around T Tauri stars. 
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3.2. 12 CO J=2-l 

Figure la shows the CO J=2-l line emission distribution over the velocity range with 
significant emission. The images show an approximately symmetric pattern around the 
central velocity of 6.8 km s _1 . At the extremes, the emission is centered at the continuum 
peak position, while at intermediate velocities the emission shows a gradient oriented slightly 
east of north that at the central velocities separates into two peaks. Figure 2 shows the 
spectrum of the CO J=2-l line convolved to a beam size of 3%, which encompasses the 
full spatial extent of the emission. The spectral profile is double peaked. All of these 
features taken together strongly suggest a disk at low inclination to the line-of-sight in 
Keplerian rotation, as calculated by Beckwith & Sargent (1993), where low J CO lines have 
high optical depth and trace gas temperature near the disk surface. A Keplerian velocity 
field is appropriate for the estimated low disk mass relative to the central star mass, i.e. 

0.02 Mq) << Mgtar (~ 2.0 Mq). 

3.2.1. Disk Model 

In order to verify the Keplerian nature of the observed rotation of the outer disk material, 
and to provide quantitative constraints on the basic disk parameters, we performed a series 
of radiative transfer calculations based on a physical model selected from the catalog of 
irradiated accretion disk models of D’Alessio et al. (2005). We selected the model that Dent 
et al. (2005b) found to closely match the spectral energy distribution of HD169142. In this 
model, the central star parameters are T e ff = 9000 K, M* = 2.0 M 0 , i?* = 1.7 R e , and age 
10 Myr, which are well matched to the stellar properties of HD169142. The disk has a mass 
accretion rate of 1.0 x 10“ 8 M 0 yr _1 , with viscosity parameter, a = 0.01. The dust in the 
disk is well mixed with the gas and has a power law size distribution n(a) ~ a~ p , p = 3.5, 
with sizes ranging from a min = 0.005 /jrn and a max = 1 mm. 

We identified an optimal disk radius, inclination, and orientation on the plane of the 
sky, and CO abundance by performing a chi-squared minimization, comparing the SMA 
visibility data and visibilities derived from a grid of disk models spanning a range for each 
parameter. The radiative transfer calculations were done with the 2D accelerated Monte 
Carlo code of Hogerheijde and van der Tak (2000) to generate simulated data cubes from 
which visibilities were derived at the ( u , v ) points observed by the SMA for the calculation 
of the chi-squared metric. The (dust) temperature and density values from the model were 
re-binned and truncated for a range of outer disk radii from 150 to 400 AU, based on the 
extent seen in the CO J=2-l images. Because of beam dilution, these data do not constrain 
the inner disk radius. The narrow line width suggests a viewing angle close to face-on for 
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the disk, and we explored a range of inclinations from 0° to 30°. The turbulent component 
of the velocity held was fixed at the low value of 0.12 km s~ 4 , less than the resolution of the 
data. The CO abundance is characterized by a depletion factor from the nominal dark cloud 
abundance of 10” 4 with respect to molecular hydrogen. 

We find a best fit disk model with radius 235 ± 5 AU and inclination 13 ± 1°, oriented 
at position angle 5 ± 5°, with a CO abundance of 2.5 x 10 -5 (depletion factor 4). We assign 
uncertainties to the fitted parameters that are somewhat larger than the formal uncertainties 
based on visual examination of the model images and their residuals. Table 2 summarizes the 
model parameters. Figure lb shows the velocity channel maps from simulated observations 
of the best fit model. The Keplerian disk model clearly captures the main features visible in 
the data. Figure 2 shows the spectrum at 3'/5 resolution derived from the model, overlayed 
on the data. This comparison gives an impression of differences integrated over the full 
spatial extent of disk emission. The derived low inclination is compatible with the lack 
of detectable ellipticity in the 1.3 millimeter dust emission, as well as the approximately 
circular appearance of the polarized scattered optical light nebula. The derived disk outer 
radius of 235 AU is larger than the 130 AU inferred from unresolved single dish CO J=3-2 
data (Dent et al. 2005a), though somewhat smaller than those found for other resolved pre- 
main-sequence A star systems, for example 545 AU for MWC480 (Simon et al. 2000) or 835 
AU for HD34282 (Pietu et al. 2003), and comparable to some lower mass T Tauri star disks 
(Simon et al. 2000). The modestly depleted CO abundance in this disk is similar to values 
typically found for the Keplerian disks around T Tauri stars, for example 1.4 x 10~ 5 derived 
for the DM Tau disk by Dutrey et al. (1997) 

The channel maps present minor asymmetries and irregularities that are not reproduced 
by the simple, smooth, Keplerian model. However, there are no dramatic departures from 
this equilibrium configuration like those seen in the circumstellar material around AB Aur, 
which is perhaps a considerably younger Herbig Ae star (Pietu et al. 2005; Corder et al. 
2005). Figure lc shows images made from the difference between the disk model and the 
data for each channel. These difference images show faint, extended emission in the data 
at blueshifted velocities, and slight excess in central peak brightness at redshifted velocities. 
These differences between the data and the model are not especially significant, but they 
may reflect real source structure. Observations with higher sensitivity and higher angular 
resolution will be needed to be more definitive. 



4. Conclusions 


We have imaged the Herbig Ae star HD169142 in 1.3 millimeter continuum emission 
and 12 CO J=2-l line emission at ~ 1'.'5 resolution. These observations resolve a disk of 
molecular gas in Keplcrian rotation around the star. The presence of the disk is strong 
evidence that HD169142 is a bona fide pre-main-sequence star, despite its distance from any 
molecular cloud. The star’s spectral energy distribution, isolation, and the regularity of its 
disk velocity held suggest an advanced pre-main-sequence evolutionary state. The millimeter 
dust emission indicates a mass of ~ 0.02 M 0 for the disk. Modeling the CO line emission 
shows that the disk has an outer radius of 235 AU and is viewed close to face-on, compatible 
with an associated scattered light nebula seen with differential polarimetry. The disk mass, 
size, and CO depletion are in line with those found for resolved Keplerian disks around other 
young stars. 


We thank Jim Moran, the SMA director, for allocating observing time to the Harvard 
University Astronomy 191 undergraduate class that allowed A.R. and M.L. to undertake 
this work. We thank Nimesh Patel for help with the observations, and Nuria Calvet for 
discussions about disk models, and for pointing us to other work underway on HD169142. 
We acknowledge NASA Origins of Solar Systems Program Grant NAG5-11777 for partial 
support. 
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Fig. 1.— (a) upper: A series of 12 CO J=2-l line velocity images of HD169142, and the 
1.3 millimeter continuum image inset at the bottom right. The contour levels for the line 
images are —2, 2,4, 6,... x 200 mJy. Negative contours are dotted. The small cross marks the 
position of the continuum peak. The contour levels for the 1.3 millimeter continuum image 
are 5,10,15,... x 5 mJy. The filled ellipse in the lower left corner of the upper left panel 
shows the 1"6 x l'.'O p.a. 20° synthesized beam, (b) middle: The 12 CO J=2-l line velocity 
images for the best fit disk model (see text). The contour levels are the same as in (a), (c) 
lower: Difference images formed from the subtraction of the best fit model visibilities from 
the observations. The contour levels are the same as in (a). 



































flux (Jy) 


HD169142 
CO J=2 —1 



Fig. 2.— Spectrum of 12 CO J=2-l line observed from HD169142 at the position of peak 
continuum emission convolved to a 3'!b beam size (histogram) and the best fit disk model 
(solid line). 
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Table 1. HD169142 Observational Parameters 


Parameter Value 


2005 April 19 (7 antennas) 
a = 18 h 24 m 19?78, S = -29 ft 46 m 49?37 
27 to 181 meters 
Callisto 


Observations 
Pointing center (J2000): 
Min/Max baseline 
Flux Calibrator 
Primary Beam HPBW 
Synthesized Beam PBW 
K/Jy 

Spectral Line setup: 
transition, frequency: 

channel spacing: 
r.m.s. (line images): 
r.m.s. (continuum) 


1"6 x 1 " 0 , p.a. 20° 

15.4 

512 channels/ 104 MHz 
12 CO J=2-l, 230.538 GHz 
0.27 km s _1 
180 mJy 
2.7 mJy 


Table 2. HD169142 Model Parameters 


Parameter 

Value 

Physical Structure 

Irradiated accretion disk a 

Velocity Field 

Keplerian Rotation 

Stellar Mass 

2.0 M 0 

Turbulent Velocity 

0.12 km s -1 

Disk Outer Radius 

235 AU 

Inclination Angle 

13° 

Disk P.A. 

5° 

CO abundance 

2.5 x 10" 5 


a D’Alessio et al. (2005) 











